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Abstract. Helioseismology has provided very detailed inferences about rotation of the solar 
interior. Within the convection zone the rotation rate roughly shares the latitudinal variation 
seen in the surface differential rotation. The transition to the nearly uniformly rotating radiative 
interior takes place in a narrow tachocline, which is likely important to the operation of the solar 
magnetic cycle. The convection-zone rotation displays zonal flows, regions of slightly more rapid 
and slow rotation, extending over much of the depth of the convection zone and converging 
towards the equator as the solar cycle progresses. In addition, there is some evidence for a 
quasi-periodic variation in rotation, with a period of around 1.3 yr, at the equator near the 
bottom of the convection zone. 
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1. Introduction 

Convection and rotation are intimately linked in the solar convection zone. The convec- 
tive timescales in the deeper parts of the convection zone are similar to the solar rotation 
period, giving rise to a strong influence of rotation on the convection dynamics. Also, 
transport of angular momentum by convection is likely important for the solar surface 
differential rotation and the interplay between convection and rotation presumably also 
causes other dynamical phenomena in the convection zone, such as the meridional circu- 
lation. Finally, interaction between convection, rotation and other flows are assumed to 
lead to the generation of the solar large-scale magnetic fields and their cyclic variation. 

Before the advent of helioseismology, observation of solar rotation was limited to the 
motion of features across the solar disk, or measurements of the Doppler velocity of the 
solar surface plasma. This showed the variation in rotation with latitude, from a rotation 
period around 25 days at the equator to around 35 days near the pole, although the 
motion at high latitudes was difficult to follow. Also, differences in rotation rate between 
the surface plasma and magnetic features presumed to be anchored at different depths 
gave some i ndication of an increase of angular v elocity with depth in the near-surface 
region fe.g.. Iwilcox fc Howa,rdllT97fl lFoukallll97l. 

Theoretical modelling through hydrod ynamical simulations of the convection zone 
l)datzmaierl Il98<4 lOilman fc Millerl Il986) was able to reproduce the surface variation 
and tended to show an angular velocity depending only on the distance from the rotation 
axis; such 'rotation on cylinders' was expe cted from simp le hydrodynamical arguments, 
leading to the Taylor-Proudman theorem l|Pedloskvlll987|) . Modelling of rotation of the 
deep solar interior was extremely uncertain. From observations of other stars it is as- 
sumed that the Sun was rotating substantially more rapidly in its early life, with a loss 
of angular momentum through the magnetized solar wind coupled to the convection zone; 
however, it was unclear to what extent transport of angular momentum in the radiative 
interior would lead to a spin-down of the solar core. 
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The observatio nal situation has been changed dramatically by helioseismology. Already 
early analyses bv lDeubner et al\ l)l979|) of high-degree five-minute oscillations showed the 
effect of the advection of the modes by rotation, providing an in dependent indication of 
an increase of the velocity with depth. iGough fc Toomrel Jl983l) provided a more solid 
theoretical basis for the extraction of such subsurface velocity signals or sound-speed 
fluctuations fro m helioseismic data. With observations over the last decade from the 
GONG network faarvev et a/.lll996j) and the MDI instrument on the SOHO spacecraft 
i Scherrer et a^J Il995^ we now have a detailed picture of rotation in the solar interior 
and its variation over a full 11-year cycle. Interestingly, this does not conform with the 
early models of rotation on cylinders in the convection zone; as discussed by Brun et al. 
(these proceedings) models of the convection zone are still not fully able to account for 
the inferred rotation profile. A striking feature is the sharply localized change from the 
latitudinally differential rotation in the convection zone to the nearly solid-body rotation 
in the radiative interi or. This takes place in a narrow region, the so-called tachocline 
I Spiegel fc Zahnlll992T ) near the base of the convection zone, of likely substantial impor- 
tance to the generation of the solar magnetic cycle. 

Here I give a brief overview of the helioseismic analyses that have led to these inferences 
of solar internal rotation and discuss some of the results. A more extensive review o f 
solar rotation and helioseismic investigations of it was given bv lThoniDSon et all <|2003), 
while IChristensen-Dalsgaard fc Thomnso discussed the observed properties of 

the tachocline region in more detail. A recent detailed review of observa t ion an d modelling 
of the dynamics of the solar convection zone was provided bv iMieschl ((2005) . 



2. Helioseismic inferences of rotation 

As a background for interpreting the helioseismic results a short introduction to he- 
lioseismolog y is probably usefu l. E xtensive treatm ents o f stellar oscillations have been 
provided bv lUnno e£ alfl|l989() and lGoughl l)l993j) . while IChristensen-Dalsgaardl l|2002h 
reviewed the techniques and results of helioseismology. 

2.1. Properties of solar oscillations 

Solar oscillations are believed to be excited stochastically by near-surface convection. 
They are observed at periods between roughly fifteen and three minutes, with largest 
amplitudes near periods of five minutes. The amplitudes per mode are minute: at most 
around 20 cms -1 in radial velocity and below a few parts per million in intensity. Even 
so, it has been possible to determine the mode frequencies with extremely high accuracy, 
thus forming the basis for the helioseismic analyses. 

A mode of solar oscillations depends on co-latitude 9 and longitude (j> as & spherical 
harmonic, characterized by the degree I and the azimuthal order m, with |m| ^ I. In 
addition, the mode is characterized by its radial order n which provides a measure of the 
number of nodes in the radial direction. The angular frequency io n i m in general depends 
on all three wave numbers. 

The observed modes of solar oscillations are acoustic (or p) modes or, at relatively high 
degree, surface gravity (or f) modes. The p modes are trapped between the photosphere 
and an inner turning point at a distance r t from the centre, determined by 

^ = , (2-1) 

where c is the adiabatic sound speed. Thus high-degree modes are trapped near the sur- 
face while low-degree modes penetrate to the solar core. For the f modes the displacement 
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Figure 1. Contour plots of rotational kernels K n i m in a solar quadrant. The modes all have 
cyclic frequencies u)/2n near 2 mHz; the following pairs of (I, m) are included: a) (5, 2); b) (20, 8); 
c) (20,17); and d) (20,20). The dotted circles mark the locations of the lower radial turning 
point r t (cf. eg. 12.11 and the dotted lines show the latitudinal turning points ipt (cf. ea. l2.2H . 

amplitude decreases exponentially with depth below the surface with an e-folding dis- 
tance of approximately Rq/1, Rq being the solar radius. The extent of the eigcnfunctions 
in latitude is determined by the properties of the spherical harmonics; asymptotically, a 
mode is confined between latitudes of ±ipt , given by 



COS Ipt 



1 + 1/2 



(2.2) 



thus modes with m ~ I are confined near the equator, whereas modes with m/l <C 1 
extend over all latitudes. 

The Sun is rotating so slowly that the centrifugal force or other higher-order effects 
of rotation can be neglected. For the observed modes the dominant effect of rotation 
on the frequencies is simply the advection of the modes by rotation: the dependence 
of an oscillation on longitude and time t can be written as cos(uj n i m t — m</>), i.e., for 
m =/= behaving as a wave running in the longitude direction; thus the frequencies of 
modes travelling in the direction of rotation are increased and frequencies of modes 
travelling in the direction opposite to rotation are decreased. It is plausible that the 
resulting frequencies are given by Lo n i m — ^>nio + rn(Q), where (f2) is a suitable average, 
determined by the properties of the mode, of the angular velocity f2(r, 9); as indicated, 

is in general a function of the distance r to the centre and 9. A more careful analysis 
shows that the rotational splitting can be written as 

Sujnim = uj nlm - uj n i — m I / K n im [r, 9)tt(r, 9)rdrd9 , (2.3) 

Jo Jo 

where the kernels K n i m are determined by the structure of the Sun and the properties 
of the eigenfunctions. It is important to note that the kernels are symmetrical around 
the equator; thus the rotational splitting is only sensitive to the similarly symmetric 
component of Q(r, 9), and analysis of rotational splittings of global modes provides no 
information about the antisymmetric component of Q. Some examples of kernels are 
illustrated in Figure ^ it is evident that they largely follow the behaviour expected from 
the turning points in r and 9. 

2.2. Inverse analysis 

By observing a large selection of modes of such varying extent it is possible to make 
combinations of the rotational splittings that provid e localized measures of f2(r, 9). Tech- 
niques for such inverse analyses were discussed bv ISchou et all l|l998jh In many cases 
the inference of fl(/'o, 9q) at some point (ro, 0o) is obtained as a linear combination of the 
data which, for example, could be taken on the form m~ 1 5(jj n i m : 

O(r o ,0o) = 'V.CrUmiro^m-^Unhn = I ° [ IC{ro,9o,r,9)n(r,9)rdrd9 , (2.4) 
nlm Jo Jo 
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Figure 2. Inferred rotation rate Q/2n in a quadrant of the Sun, obtained by means of inversion 
of 144 days of MDI data. The equator is at the horizontal axis and the pole is at the vertical 
axis, both axes being labelled by fractional radius. Some contours are labelled in nHz, and, for 
clarity, selected contours are shown as bold. The dashed circle is at the base of the convection 
zone and the tick marks at the edge of the outer circle are at latitudes 15°, 30°, 45°, 60°, and 
75°. The shaded area indicates the region in the Sun where no reliable inference can be made 
with the prese nt data. The slant ed dotted lines are at an angle of 27° with the rotation axis. 
(Adapted from lSchou et qZ.lll993) . 

by using equation (|2.3[l . Here the averaging kernel is given by 

K.(r o ,6 o ,r,0) = J~] c»; ro (r , 0o)K n i m (r, 6) ; (2.5) 

nlm 

it is typically denned such as to have unit integral over (r, 8) . Thus equation (|2.4|l defines 
the inferred angular velocity as an average which ideally is localized near (ro,0o); the 
degree of localization is determined by the extent of the averaging kernel. Also, the error 
in O(ro,0o) can obviously be determined from equation l|2.4|l if the error properties of 
the data are known. How the inversion coefficients c n i m (ro, 9q) are obtained depends on 
the properties of the specific inversion method; in some techniques they are found so as 
explicitly to localize the averaging kernel, whereas other techniques attem pt to fit the 
inferre d angular velocity to the observations in a least-squares sense (see ISchou et all 
l|l998l) for details). 



3. The solar internal rotation 

The inversion techniques discussed above hav e been extensively a pplied to observed so- 
lar oscillation frequencies. Already early results ijDuvall et aZjl98^) clearly demonstrated 
that the radiative interior rotated at roughly the surface rate, with little or no indication 
of a rapidly rotating core. With the determination of rotational splittings for different val- 
ues of m it became possible to o btain information about the variation of the latitudinally 
differential rotation with depth ( Brown fc Morrowli987H IChristensen-Dalseaard fc Schoul 
1988; Br own et a/.lll989tlDziembowski et alJll989j) . demonstrating the transition between 
differential rotation in the con vection zone a nd ne ar-constant rotation in the radiative 
interior. Extensive analyses by ISchou et all l|l998j) have showed a considerable degree 
of consistency between results of different inversion methods, although some sensitivity 
remains to the techniques used in the basic analysis to determine the rotational splittings 
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Figure 3. Averaging kernels for the inversion shown in Figure^ targeted at the following radii 
and latitudes in the Sun: (0.55i?©, 60°), (0.7R@, 0°), (0.7R e , 6 0°), and (0.95Jfo 60°). The 
corresponding locations are indicated with crosses. (Adapted from lSchou et aZ.lll99Sft . 

from the observations l|Schou et al\\200^ . particularly at high latitudes. Using also data 
on low-degree modes from unresolved observations the absence of a rapidly rotating core 
has been confirmed, althou gh the inferences ar e still uncertain in the inner around 20 
% of the solar radius (e.g.. IChaolin et flHll999l) . This confirms that an efficient mecha- 
nism must have linked the radiative interior to the loss of angular momentum from the 
convection zone; the nature of this mechanism is still debated, however. 

3.1. The bulk of the convection zone 

Here I concentrate on rotation within and just below the convection zone. Figure |21 shows 
an inferred rotation profile in the outer parts of the Sun. Some averaging kernels for this 
inversion are illustrated in FigureEl evidently the combinations of the observed splittings 
are able to determine the angular velocity with fairly high resolution, particularly in the 
radial direction near the surface. In the bulk of the convection zone the latitudinal varia- 
tion clearly roughly follows the variation at the surface, the helioseismically determined 
values close to the surface be ing essentia lly consistent with direct Doppler measurements 
of the surface rotation fe.g.. Ir3eckll2000l) . It is obvious that the rotation is not constant 
on cylinders, as simple models predicted. As noted bv lGilman fc Howel l)2003l) it rather 
appears that in much of the convection zone rotation is constant on lines inclined at 27° 
to the rotation axis; the physical reason for this property is not understood. 

Very near the surface the contours clearly indicate that the rotation rate decreases 
outwards, at low and intermediate latitude. This corresponds to the behaviour inferred 
from surface marke rs presumed to be anc hored at different depths and is reminiscent of 
the early results oflDeubner et al\ l|l979l) . In a more careful analysis based on f modes 
ICorbard fc Thompson! ((20021) determined the gradient of the angular velocity in this 
near-surface region, as shown in Figure 0] Interestingly, dlnf2/dlnr ~ —1 at latitudes 
below around 40°; this is not consistent with the, perhaps natural, assumption of constant 
angular momentum maintained in this region by convective motions, which would instead 
correspond to dlnil/dlnr = —2. 

3.2. Properties of the tachocline 

As already mentioned, the properties of the tachocline are of very considerable interest. 
Particularly important are the location and thickness of the transition, first addressed 
in detail by iKosovichevI l^)9<5) . Determination of these properties must obviously take 
into account the resolution of the inversion, as characterized by the averaging kernels 
(cf. Figure 0J. The effect on the inversion is illustrated in Figure^, on the basis of 
inversion of artificial data computed for the rotation profile shown by the heavy curves. 
The corresponding inversion results, indicated by thin curves and shaded areas, clearly 
show a transition that is more gradual than for the 'true' underlying rotation rate, 
as a result of the smearing by the averaging kernel. The result of the corresponding 
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Figure 4. Logarithmic gradient of the angular velocity obtained from linear fits in the outer 
two per cent of the solar ra dius to the results of inversion of f-mode splittings. (Adapted from 
iCorbard fc Thompsonll2003) . 



inversion of solar data l|Charbonneau et oi.lll99ffl . shown in Figure 03, in fact indicates 
a width that roughly corresponds to the width assumed for the artificial data. Further 
information about the transition can be obtained, as discussed by Charbonneau et ai, by 
arranging the linear combination of the data in equation (|2.4|l such that the corresponding 
combination K, of the kernels provides an average of the radial gradient of Q. The result 
(Figure Et) clearly shows the localized gradient in the tachocline, furthermore hinting 
that the transition takes place closer to the surface at latitude 60° than at the equator. 

From tests on artificial data, and fro m parametrized fits, it is to some extent possibl e 
to correct for the finite resolution (e.g.. ICharbonneau et aZJIl999t ICorbard et ctZJIl999l> . 
and hence obtain a measure of the actual width. For example, iBasu fc Antial l|2003|) 
determined the location r c of the midpoint of the transition and its width u>f as r c = 
O.692i? , w = 0.033i? Q at the equator and r c = 0.710i? o , w = O.O76i? at latitude 60°. 
Thus they confirmed the prolate nature of the transition and furthermore found that 
it is broader at high latitude than at the equator. For comparison, the radius at the 
base of the convection zone, assu ming spherical symmetry, has been determined from 
helioseismology as 0.713i?© (e.g.. IChristensen-Dalsgaard et oi.lll99ll) : thus, under this 
assumption, most of the transition takes place below the convection zone at the equator, 
whereas at higher latitude there is considerable overlap between the tachocline and the 
convection zone. 



4. Variations with time 

Doppler observations of the solar surface rotation have show n variations with time in 
the so-called torsional oscillations l|Howard fc LaBontelll980]) . bands of slightly faster 
and slower rotation that converge towards the equator as the solar cycle progresses. The 
availability of detailed helioseismic observations extending over a sunspot cycle now al- 
lows investigations of variations in solar structure, rotation and other types of dynamical 
phenomena beneath the solar surface. Particularly detailed investigations have been made 
of changes in rotation; the principle is to analyse the data in segments of typically 2-3 
months, and consider the residual from the time-averaged rotation rate , as a function of 
r and 9. A recent overview of the results was provided bv lHowelll2n0rt . 



f defined such that 84 % of the transition takes place between r c — w and r c + w 
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Figure 5. (a) Inversion of synthetic dataset for artificial rotation profile, based on artificial 
data with properties corresponding to an observed dataset. The thick curves show the underlying 
input rotation profile, in the form of radial cuts at latitudes of 60° (dashed curve), 30° (dotted 
curve), and 0° (equator; solid curve). The thin curves show a typical inversion solution, with 
the la range, determined from the errors in the observed mode set, indicated by the shaded 
areas, (b) Corresponding inversion result for the solar data, (c) Results of inversion for the radial 
gradient of the solar rotation rate at latitude 0° (circles), 30° (asterisks) and 60° (diamonds); 
here the vertical error bars indicate la errors, whereas the horizontal bars indicate the resolution 
of the inversion as determ ined by the extent in radius of the averaging kernels. (Adapted from 
ICharbonneau et aZ.lll999f) . 

4.1. Tachocline oscillations? 

Perhaps the most surprising result of such analyses was the detection of what appeared 
to be an oscillation with a period of 1.3 yr in the equatorial rotation rate near the base 
of the convection zone ijHowe et al\ l2000a|) . As shown in Figure H3 this oscillation was 
present, with largely consistent behaviour, in analyses of two different datasets and using 
two different inversion methods, between late 1995 and 2000, while the variation has been 
more erratic since then. The early period also appeared to show an oscillation with a 
similar period, but the opposite phase, at somewhat greater depth at the equator. Given 
the subtle nature of the signal, it is important to check whether it could be an artefact 
of, for example, va riations in the select ion of mod es included in the ana lysis; this appears 
not to be the case l|Toomre et al 1l2003|) . However. lBasu k Antial l|200l|) . while apparently 
detecting a similar variation, questioned its statistical significance. It is evident that the 
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Figure 6. Residuals from average rotation rate at r — 0.72Rq at the equator. Triangles and 
circles are based on MDI and GONG data, respectivel y, and open and c losed symbols correspond 
to two different inversion techniques. (Adapted from Howe 



erratic nature of the variation may cast doubt on its physical reality; on the other hand, 
it is not implausible that such behaviour, of possibly magnetic origin, could be linked to 
specific phases in the solar cycle. Thus continued observations during the rising phase of 
the next cycle will be of great interest. 

4.2. Zonal flows 

Helioseismic investigations have showed that, far from being a supe rficial phenomenon, 
the torsional oscillations seen on t he solar surface ifrloward & LaBontell980l) extend quite 
deeply. From inversion of f modes iKosovichev fc SchovJ l|l997l) found that the regions of 
slightly faster an d slow e r rota tion extended a few per cent of the solar radius beneath 
the solar surface. ISchonl (1999) considered f-mode data covering more than two years and 
found that these subsurface zonal flows shared the propagation towards the equator seen 
on the surface, closely linked to the 'butterfly dia gram' of the sunspots . Mo re extensive 
data, o n both p and f modes, were considered bv lAntia fc BasvJ (|200Cf) and lHowe et all 
who confirmed the equator-ward propagation of the zonal flows at low and 
intermediate latitude, the flow extending through up to one third of the convection- 
zone depth. At high latitude, on the other hand, flows of somewhat higher amplitude 
propagated towards the pole . Even deeper penetration, particularly at high latitude, was 
inferred bv IVorontsov et all l|2002h . using a non-linear inversion technique. They also 
found that the propagation could be fitted as having an 11-yr period and determined the 
amplitude and phase of this variation. 

As an example of these results, Figure compares flows from Doppler observations 
at the solar surface with helioseismically inferred flows at a depth of O.Oli?©. The two 
patterns are evidently very similar. The variation with depth is illustrated in Figure |SJ 
this confirms that the flows can be followed through most of the convection zone. Also, 
strikingly, the flow pattern appears to propagate towards th e solar surface as time pro- 
gresses; this is difficult to reconcile with a model proposed bv lSpruitl 1 200.^ according to 
which the flow originates from thermal effects at the solar surface. 

The origin of t hese zonal flows has been discussed in the context of mean-field dy- 
namo models (e-g.- ICovas et aJl200(ilRempe]|2006albli . Rempel found that the pole- ward 
propagating branch could be explained by the feed-back of the Lorentz force on differ- 
ential rotation in such models. He also showed that a similar mechanical model for the 
equator-ward propagating low-latitude branch would be strongly affected by the Taylor- 
Proudman theorem, leading to a variation with depth and latitude that is inconsistent 
with observations; thus he concluded that thermal effects must be involved. 



Rotation of the Solar Convection Zone 



127 




1995.0 1997.2 1999.4 2001.6 2003.8 2006.0 
Date (years) 




Figure 7. Zonal flow patterns, obtained as residuals from time-averaged rotation rates. The 
top panel shows the 'torsional oscillations' obtained from Mount Wilson surface Doppler ob- 
servations. The two lower panels are from inverse analyses targeted at r = 0.99-Rq; in the left 
panel global inversion of MDI data was used, whereas the right panel is based on a local he- 
lioseismic analysis with the ring-diagram technique (see Section during time periods where 
high- resolution MDI data were available. (Adapted from iHowe e~a*j l2006a 
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Figure 8. Rotation residuals as a function of time and fractional distance from the centre, 
from inversions of MDI data, at latitudes (left to right) 0°, 15°, 30°, 45° and 60°. The data have 
been averaged over periods of 1 yr. Darker regions c orrespond to rotatio n slower than average, 
and lighter regions to faster rotation. (Adapted from [Howe_e£ al. 2006b) 



5. Local helioseismology 

I have so far considered global helioseismology, based on frequencies of global modes 
of solar oscillation. This is restricted to consider only rotationally symmetric aspects 
of the solar interior structure and dynamics, such as rotation; furthermore, as noted in 
Section 12.11 the analysis is sensitive only to aspects that are symmetrical around the 
equator. These restrictions can be avoided through the use of local helioseismology where 
the data are analysed instead in terms of wave propagation in a smaller part of the solar 
surface. An extensive review of the techniques and results of such analyses was given by 
iGizon k Bird] l|2005l) . 

Detailed investigations of sub-surface flows have been made with the ring- diagram 
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technique, following the early analysis of iGough k. Toomrel l|l983t) . The oscillations are 
analysed in smaller patches on the solar s urface, to produce local power spectra in terms of 
the horizontal wave vector k h llHil]lll988l) . By considering a large number of such patches 
distributed over the solar disk one can build up a picture of the subsurface structure and 
flows as a fun ction of position o n the Sun. A detailed investigation of this nature was 
carried out by Ha ber et they determined the longitudinally averaged rotational 

flow, also identifying the zonal flows. The results were similar to those of the global 
analyses, but with clear differences between the northern and southern hemispheres. 
This is also illustrated in Figured where the lower right p anel shows results obtained 
with a ring-diagram analysis. Similar results were obtained bv lZhao fc Kosovichevl l|2004) 
using the time- distance technique, where a correlation analysis is used to infer the travel 
time of waves between different points on the solar surface. 

Local helioseismology has also revealed other kinds of large-scale flows beneath the 
solar surface. An important example is the meridional cir culation, which is also seen 
on the solar surface in Doppler-velocity observations (e.g.. lHathawavlll99fil), such that 
the flo w is generally in the pole- ward direction. Using time-distance analysis Idles et "all 
(1997) found a si milar flow in t he up per parts of the convection zone . This was analysed 
in more detail bv lHaber et al\ (|002) and lZhao fc Kosovichevl (|2004^> using ring-diagram 
and time-distance analyses, respectively. On smaller scales local analyses have revealed 
convective flows on supergranular scales, as well as complex flows associated with active 
regions (see also Kosovichev, these proceedings). 



6. Open questions 

While major progress has obviously been made on our knowledge about rotation in 
the solar interior, much remains to be done. The determination of the rotation of the 
inner solar core is still highly uncertain. Unfortunately, substantial improvement in the 
low-degree p-mode splittings, required to determine the core rotation, will require much 
more extensive observations than now available; observation of g modes, which are more 
sensitive to the properties o f the solar core, remains elusive despite recent encouraging 
progress ((Garcia et a/.ll20f)fil) . Better observational constraints are also required on the 
properties of the tachocline region; this includes details of solar structure such as the 
depth of the convective envelope and the detailed properties of the transition to the 
radiative region, as well as a better determination of the location and thickness of the 
transition in rotation. Also, it is obviously important to test the reality of the 1.3 yr 
oscillation in rotation in the tachocline region and investigate its diagnostic potential. 

From a theoretical point of view, the evolution of rotation in the radiative inte- 
rior of the Sun is still poorly understood, although there are models that are able 
to reproduce the present near-co nstant rotation through , for example, the effect of a 
dynamo-generated magnetic field ( Eg genberger et fli]l2005l) or transport by gravity waves 
( Charbonnel & Talonll20f)5^ . This issue is obviously closely related to the establishment 
of the tach ocline, where magnetic fie lds have also been invoked as the most likely mecha- 
nism (e.g.. lGough fe Mclntvrdll9 98). Progress towards understanding the mean rotation 
profile within the convection zone has been ma de in terms of mean-field models which 
involve thermal effects in the tachocline region l|Rempe]|l2005|) . Evidently, the preferred 
route to understand the dyna mics of the conv e ction zone would be through realistic 
three-dimensional simulations. iBrun fc Toomrel ((2002) showed that increased numeri- 
cal resolution in such simulations improved the agreement with the observed profile, 
and the simulations have confirmed that including entropy variations in the tachocline 
region can act to bring the results into better agreement with the helioseismic infer- 
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ences l|Miesch et alJEooil see also Brun et al. th ese proceedings). Int erestingly, detailed 
modelling of the outer 5 % of the solar radius ()DeRosa et aZ.ll2002^ has reproduced a 
near-surface decrease in the angular velocity, such as observed in the shear layer (cf. 
Figure A definite model of the zonal flows is still lacking; one may hope that such 
flows will eventually emerge as natural features of the hydrodynamical simulations. 
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Discussion 

Pecker: The trend towards the equator seems to be much quicker (in your diagram) 
for torsional oscillations than it is for helioseismological data. Why is that so? 

Christensen-Dalsgaard: In fact, my impression from Figure0is that the helioseis- 
mically inferred flow is quite similar to the behaviour of the surface torsional oscillation 
and, if anything, is at slightly lower latitude. This could be related to the apparent 
propagation towards the surface, as illustrated in Figure |H1 



